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Abstract 

The lepton flavor violating process J/ip — » II' (I ^ V) serves as an ideal place to probe the unparticle 
theory. Such process can only occur at loop level in the Standard model (SM), so that should be 
very suppressed, by contrast in unparticle scenario, it happens at tree level and its contribution 
may be sizable for practical measurement. Moreover, the BESIII will offer the largest database 
on J ftp which makes more accurate measurements possible. Furthermore, for such purely leptonic 
decays background is relatively low and signal would be cleaner. Our work carefully investigates 
the possibility of observing such processes from both theoretical and experimental aspects. 
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I. INTRODUCTION 



Violation of the lepton flavor symmetry plays an important role for exploring new physics 
beyond the standard model (SM). Neutrino oscillations which are observed in many experi- 
ments, indicate that lepton flavor symmetry is broken and the minimal SM must be extended 
to accommodate non-zero neutrino masses. Up to present, the lepton flavor violation (LFV) 
has only been observed in the neutrino sector. This violation is induced by the neutrino 
masses and existence of the right-handed neutrinos. By contrast, for the charged lepton 
sector, there is no any evidence for such a violation. Many studies have been devoted to 
explore the LFV through r and Z decays In Ref. (2|, it was suggested to look for 
LFV via charmonium J /if} decays. In terms of the database of 5.8 x 10 7 J /if? data at BES 
processes J /if) — > e// J /if) — > lit, er [5j have been explored and upper limits for their 
branching fractions were set. New hope is raised that the forthcoming upgraded BESIII 
will accumulate up to 10 10 J /if) per year [fj], which makes it possible to test LFV with a 
much higher precision. It is well known that the SM contribution to such violation can only 
occur via loops and are very suppressed, so that cannot produce observable effects even with 
so large database of J/if>. Therefore it is time to investigate possibility of discovering such 
violation in terms of the theoretical models beyond SM. 

Recently, Georgi proposed an interesting idea that a scale invariant stuff may exist in 
our world. The scale invariant stuff contains no particle, but the so-called unparticle [3]. 
To be consistent with the present experimental observations, in the effective theory about 
the unparticle, the coupling of unparticle to the ordinary Standard Model (SM) matter 
must be sufficiently weak. The scale dimension of unparticle is in general fractional rather 
than an integral number. The interactions between the unparticle and the SM particles 
in low energy effective theory can lead to various interesting phenomenology. There have 
been many phenomenological and theoretical explorations related to unparticles in literature. 
Since interaction between unparticle and ordinary SM matter is rather weak, we can expect 
that its participation cannot produce observable effects for the processes (production and 
decays) which have large rates, but definitely may play important roles in the rare processes 
where the SM contribution are suppressed. 

In unparticle physics, the unparticle can couple to different flavors of the SM leptons and 
thus induce LFV even at tree level. There have been several studies of LFV in unparticle 
physics, for instance, Jj— » 3e in js| ; ll — > e+U, 3e in jsj; M° — > IV, e + e~ — ► IV in 10]; ll — > ej, 



ll — e conversion in JJJ; muonium and antimuonium oscillation in 12]; r — > ll', I — > Vjj in 
131 ] . In this work, we propose to study LFV in unparticle physics via J /if) — > IV decays. There 
are several advantages of exploring LFV in J /if) — > IV decays, in particular the J /if) — > eii 
process. On the experimental side, the final states contain only leptons which are easy to 
detect. Compared to e + e~ — > IV processes which occur due to non-resonance contribution, a 
large database of J /if) is expected at BESIII which will begin running very soon. From the 



2 



theoretical point of view, the processes J/ip — > IV are free of hadronic uncertainties after the 
parameter related to the binding effect of J/ip is experimentally determined by measuring 
J/ij) — > e + e~({i + fi~). 

Anyhow, an observation of J/ip — > IV (I 7^ /') would be a clear indication of new physics 
beyond SM. In our case, cc in J/ip annihilate into a virtual unparticle which later converts 
into two leptons IV with different flavors. As we will show that the theoretical prediction on 
the branching ratio of J/ip — > e/i can reach about 10~ 8 order, which is absolutely possible 
to be observed at BESIII. To help or experimental colleagues, on the basis of BESIII envi- 
ronments, we carry out an investigation on the possibility of search for J/ip — > e/i in terms 
of a Monta Carlo simulation. 



II. J/Y> -> IV IN UNPARTICLE PHYSICS 

We start with a brief review about the unparticle physics which is relevant to this study. 
The scale invariant unparticle fields emerge below an energy scale A u which is assumed to be 
at the order of TeV. The interactions of the unparticle with the SM particles are described 
by a low energy effective theory. The coupling of unparticle to SM fermions (quarks and 
leptons) is generally given by the following effective operators as 

4hf'fOu, ^r-rtfOu, %f\fd»o u , §;fi^o u - 

ly U ly U ly U Ly u 

#mo£, jgr-rsfdM, -^r/Vo& t#i/W^, (i) 

iy u iy u iy u iy u 

where Ou and O u denote the scalar and vector unparticle fields, respectively. Tensor un- 
particle does not contribute in our case. The are dimensionless coefficients. Because J/ip 
is a vector, only the vector current f'j^f couples to J/ip. For a scalar unparticle coupling 
f'lufd^Oui it is proportional to the momentum of the unparticle. Using the equations of 
motion, one can immediately show that this part is proportional to the lepton masses in 

decays of J/tp — > IV and thus is suppressed. Consequently, we will only consider the effective 

/'/ 

interaction ^-1 J^JO^ which dominates in the decays of J/ip — > IV. 
For the vector unparticle field, the propagator is given by 

A n iiv 1 pfi pv I p2 

^<0TOx)P&(0)|0> = ^ m ^\ p : +tey i e-^. (2) 

16tt 5 / 2 T(d u + 1/2) 
du (2ir)^T(d u -l)T(2duY {) 

The function sin(<i^7r) at the denominator implies that the scale dimension du cannot be inte- 
gers except du — 1. The peculiar unparticle propagator and the phase factor e~ l ^ du ~ 2 ^ n which 
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where 
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provides a CP conserving phase produces novel effects in high energy scattering processes 
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CP violation 15j, neutral meson mixing [16| and neutrino processes 17] 

/ 




FIG. 1: The diagram for the decay of J/ip — > IV . The double dashed lines represent the unparticle. 

The lowest order contribution to decays of J/ip — > IV is a simple tree diagram, which is 
depicted in Fig. [IJ From the diagram, the exchange of vector unparticle is analogous to 
a photon exchange unless photon cannot couple to flavor changing currents. By using the 
effective operators and unparticle propagator, it is straightforward to write out the decay 
amplitude as 

c vv c vv A d u m ^V ? -,f/ r y, v (u\ 
A^- 2 2sin^vr W" 1 * h 

LA 



M {J/ip{P) -> l{k)V{k')) 



(4) 



where s = mi. The decay constant is derived from the pure-leptonic process J/ip — ► 
e + e~(/i + /i~). In principle, unparticle can also contribute to J/ip — > e + e~, but compared 
with the QED contribution, its fraction in BR(J/ip — > e + e~) ~ 6% and BR(J/ip — > e/i) < 
1.1 X 10~ 6 3] is completely negligible. 

The decay width of J/ip — ► ZZ' is given by 



r(j/^ -> z+z'-) 



spins 



with 



1^ ^iCC s*it& 
I c vv c vv 



spins 



A. 



du 



3 V A^" 



2smdu7i s 



1-d u 



2(P • k){P ■ k') + Jfe • k' 



(5) 



(6) 



For J/ip — > e/i, where the final lepton masses are both small compared to mass of J/ip, 
we have the simplified form for the decay ratio as 

j ,2 



r(J/V> -> e+/i~) 
r(J/^ -»• e+e") 



47rae c 2sincZ^7r V^-w 



(7) 



with e c = 2/3. Noting that T(J/ip 
have 



e /i + ) = T(J/ip — > e + /i ) due to CP symmetry, we 



T(J/ip -> e/i) = r(J/V> -> e-/i+) + r(J/V> -> e+/i") 
= 2r(J/^ -f e + /0- 
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We would like to present several comments as follows: (1) The relation of Eq. ([7]) can be 
applied to other vector meson whose mass is larger than that of J/ip, such as T, etc.. (2) 
The decay width of V — > efi depends on the invariant mass squared s, and is proportional 
to (s/A^) dw_1 . When du > 1, the decay width decreases as a power function. (3) For 
1 < dy < 2, the decay width decreases faster as du increases. The dependence of the 
branching ratio of J/ip — > efi) on dimension du is shown in Fig. [2J 

An important constraint on the unparticle interaction which is related to J/ip — > IV pro- 
cesses comes from the invisible decay of J/ip. In SM, the invisible final states are neutrinos. 
The J/ip — > vv decays are highly suppressed and can be neglected. In unparticle physics, 
the unparticle U cannot be detected by experimental apparatus and if J/ijj decays into an 
unparticle pair, Ws would correspond to missing energy. Moreover, the J/ip — > U decay 
is possible because the mass of U is not fixed, and IA also escapes from our detector. It 
provides a constraint on the the coupling of unparticle to cc quark pair. The rate of the 
J/ijj — > Li decay are derived using the same effective interactions as that in J/ip — > ll f , and 
it is straightforwardly to obtain 

T{J/1)^U) _ 3A du \c c v c v \ 2 ( s 
r(J/<0 -> e+e-(/x+/r-)) 8na 2 e 2 c \Al 

Our result is in agreement with [18| except a factor of 4 difference due to different conventions. 
Another constraint on the coupling of unparticle to e/i may come from fi —>■ e+U and fi —> 3e. 
In (sl, the scale dimension for vector unparticle d u > 2 is used. We will follow 0] and assume 
1 < du < 2. If du > 2 as in Q], the branching ratio of J/ip — > e/i will be so small that it 
cannot be observed at BESIII only if du is close to 2 (other integers are also possible). 
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FIG. 2: The branching ratio for J /ip — > e/i with scale dimension du- The dashed line represents 
the present experimental upper bound. 

The BESII collaboration recently sets an upper limit for the invisible decay ratio 



< 1.2 x 10" 



20| | . It determines the upper bound on the coupling coefficient of 



unparticle to cc quarks. If choosing Au = 1 TeV and du = 1.35, then we obtain \cy V \ < 0.01. 
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The present experimental upper bounds on the branching fractions are given in |4j, |5j: 
BR(J/?p -> e/i) < 1.1 x 10" 6 , BR(J/ip -> er) < 2.0 x KT 6 , BR(J/4> -> /xr) < 8.3 x 10" 6 . 
Among them, the upper bound on J/if> — > e/x is the lowest which may provide the most 
stringent constraint on the coupling. If we adopt the constraints from invisible decay to 
fix the unparticle parameters, i.e. du = 1.35, Au = 1 TeV and choose Cy V Cy V = 0.01, the 
branching fraction is predicted to be BR(J/ip — > e/i) = 7 x 10~ 8 . We may further constrain 
the unparticle parameters by the upper bound on BR(J/ip — > e/t). When A u and Cy V Cy V 
are fixed, for instance, choosing Au = 1 TeV and Cy V Cy V = 0.01, we obtain du > 1-26. 
This result is consistent with the constraint achieved from D° — D° mixing within la error 
tolerance 12 ]: cf^ = (0.21 + n) ± 0.07 where n is an arbitrary integer (note that uncertainty 
is still large). We will discuss the possibility of observing fractions of order 10~ 8 at BESIII 
in the next section. 

In {2)], the authors obtain a bound BR(J/ip — > e/i) < 8.3 x 10~ 6 based on a model- 
independent approach which introduces an effective four-fermion contact interaction. By 
considering the strong bounds from three body LFV processes, such as /i — > 3e and using 
the "unitarity inspired" relation between two and three-body LFV decays, it is found that the 
upper bound is constrained to be very small BR(J/if) — > e/i) < 4 x 10~ 13 [3]. However, it is 
noted that such serious bound can be avoided due to some possible kinematical suppression or 
cancelations |2]. Our predicted value lies a range between [2J and s|. The precise prediction 
is difficult because of the unknown model parameters, especially the scale dimension l . 



III. ESTIMATING POSSIBILITY OF THE DETERMINATION OF BR(J/ip -» e/i) 

BESIII is a spectrometer which will begin operation in 2008 at the upgraded Beijing 
Electron-Positron Collider (BEPCII) whose designed luminosity is approximately 100 times 
higher than that of BEPC [6| . With this high luminosity, the BESIII detector will be able to 
collect over 10 billion J /if} events in one year of running. The performance of the BESIII 



will be much more effective in comparison with the BESII 21(, so that we can expect that 
BESIII would perform a more precise and clear particle identification (PID) than BESII. It 
helps to accurately analyze rare decay channels, such as the concerned lepton-flavor violating 
processes. 

The upper bound on BR(J/ip — > e/i) was achieved with the the BESII data [4] a long while 
ago. Before BESIII begins running, we would roughly pre-estimate the upper bound with 



1 For scale dimension of the vector unparticle, it is pointed out that conformal symmetry puts a lower bound 
d > 3 in 19(. If we take this constraint, the fraction of BR(J/ip — > IV) will be decreased by more than 
two orders and it is impossible to observe them in experiment except du is very close to 3. One may avoid 
this problem by taking a point of view that scale symmetry is different from conformal symmetry. The 
discussion of this controversial topic is beyond this study. 
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TABLE I: Monte Carlo geometric efficiency. 



channel 






MC Efficiency 


J/ip —> e[i 




£ ef i-MC 


53.7% 


J/ip — > ee 




f-ee-MC 


61.47% 


J/ip -> /U/i 






58.32% 


J/tp — » 7T7T 






52.74% 


J/V> -» ^T-ftT 




tKK-MC 


24.38% 


e + e~ — > e + e~ 


(7) 


e ee(-y)-MC 


32.51% 


e + e~ — > 


"(7) 




42.96% 



the BESIII data, if signals are still not observed. For BESIII, ten billion J/ip events can be 
used for our analysis. The main backgrounds to our signal channel J/ip —>■ e/x are the nearly 
back-to-back tracks in J/ip — > e + e~, K + K~, 7r + 7r _ and e + e~ — > e + e~(7), fi + fi~~(-y). 

Since the BESIII experiment has not begun running yet, some detailed information about 
BESIII environment is not available, therefore we adopt relevant parameters for BESII to 
estimate BR(J/ip —>■ e/i) and the corresponding backgrounds, except the number of available 
J/ip events. 

The method of determining the selection efficiency for J/ip — > e\i and the rate of misiden- 
tifying the background as a signal is the same as the one used in [4|. The overall selection 
efficiency for the signal channel (or rate of misidentifying the background channel as the sig- 
nal channel) includes the /1 and e particle identification (PID) efficiency and the geometric 
efficiency. Indeed, the geometric efficiency of BESIII is only slightly different from BESII, 
thus in this work we set the same geometric efficiency as that for the previous BESII exper- 
iments. The PID efficiency and contamination rate which refer to identifying other particles 
as electrons or muons in BESIII are taken from In order to discriminate muons to reject 
other particles as much as possible, a very rigorous criterion of the /1 identification has to be 
performed while selecting signal event J/ip — > e/i. The geometric and the PID efficiency are 
listed in Tables fl] and [III respectively. 

TABLE II: The particle identification/misidentification efficiency. 





identified as e 


identified as \x 


e 


98.0% 




A* 




19.0% 


7T 


0.1% 


0.46% 


K 


0.1% 


0.38% 
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The overall selection efficiency for the signal J /if) — > e/j, (or rate of misidentifying the 
background as the signal J /if) — > e\i) can be calculated with the PID (or misidentification) 
and geometric efficiency given in the tables. For the signal channel J/if> —>■ efi, the total 
selection efficiency is then 

where e e(I _Mc is the geometric efficiency, e e ^ e is the electron PID efficiency, and e^-^ is the 
muon PID efficiency. With the number given in Tables [J and [Til we obtain the selection 
efficiency for J /if) — > e\x as about 10%. One can evaluate the misidentification rate of 
the background channels J/if> — > XX or e + e~ — > XX, where X = e, /i, tt, X, (i.e. the 
background channels are misidentified as the signal channel) as 

exx = exx-MC x ex^e x e x ^ M x 2. (11) 

where exx-uc is the Monte-Carlo geometric efficiency, shown in Table [IJ and ex^ e and 
£x-*/x are the contamination rate for X being identified as an electron or a muon. 

Because electrons and muons cannot be misidentified with each other, the only back- 
ground which should be taken into account in the estimation is the hadronic channels. The 
misidentification rate and the number of background from hadronic channels are estimated 
in terms of the methods given in [4], and the results are listed in Table [TTT1 



TABLE III: The misidentification rate and the number of background from hadronic channel. 



decays 


misidentification 


number of 




rate 


background 


J /if) — > -kit 


4.85 x 1(T 6 


7.13 


J/V> -> KK 


1.85 x 1(T 6 


4.38 


total 




11.5 



Now, we can set an upper limit for BR(J/ip e/x) if there only few events are observed 
in BESIII and they are consistent with the background estimation. For an illustration, we 
discuss a case similar to that we did for BESII where there were four J /if) — > efi candidates 
observed by the detector Although our real case may be different, the analysis should 
provide a useful reference. With the four observed events which are supposed to be the signal 
and 11.5 background events from hadronic channels in BESIII, we gave an upper limit for 
BR(J/if) — > e/i) as 

BR < X(N OB ,N BG )/[N T x ej/^ efl }. (12) 
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where A is calculated using the method of 23| to be A=2.31 is obtained at 90% C.L.. There- 
fore, we obtain BR(J/ty — > e/i) < 2.3 x 10~ 9 in the BESIII experiment, which is more precise 
than the BESII by two orders. 

Now let us discuss the case with branching ratio being about BR(J/ip — > efi) = 7 x 1CT 8 
as the unparticle physics predicts. If the total selection efficiency for J/ip — ► efi is 10%, 
then, about 70 events are passed through the selection criterion of J/tj) — > efi in 10 billion 
J/ijj events collected by BESIII during one year of running. From Table 3, there are 11.5 
background events as estimated. Then, the signal to noise ratio reaches about 7:1. It is 
possible to determinate the branching ratio of J/ip — > e/i since the signal significance is 
enough. 

In fact, the electron and muon identification is very important in the analysis of the 
concerned channels. The muon identification in BESIII is better than the one in BESII. 
Certainly, the contamination rate for kaons being misidentified as muons should be lower 
than the one for pions j2l|. In the work, however, the contamination rate for kaons be- 
ing misidentified as muons has to be set to the same as the contamination rate for pions 
being misidentified as muons, because of lack of precise information about kaons. So the 
background for J/ip — > e/i is overestimated in this work. In fact, the background in BESIII 
should be less than that estimated in the work. Therefore, the real signal to noise ratio for 
J/ip — > e/i is higher in real BESIII experiments. 



IV. CONCLUSIONS 



The lepton flavor violating decays of J/ip — > IV provide an ideal probe to explore new 
physics beyond the minimal standard model. In this study, we have presented a detailed 
analysis on the contribution of unparticle physics to J/tp — > 11' decays. The coupling of 
unparticle to different flavors of leptons cause the LFV processes at tree level. 

Even though the scaling factor (m^,/A^) dw_1 which is sensitive to the scale dimension 
of the unparticle field suppresses LFV, one still hopes to observe such rare processes at 
experiments with very high luminosity and precision. Within a reasonably chosen parameter 
space, the branching fraction of J/ip — > efi can reach about an order of 10~ 8 which is possible 
to be observed at the forthcoming BESIII. 

Indeed, the BESIII is planning to produce 10 10 J/ip per year. These huge data samples 
guarantee an exploration of LFV with a high precision. As a moderate estimation, the upper 
limit of measuring BR(J/ip — > efi) is obtained as 2 x 10~ 9 at BESIII. If so, the unparticle 
theory indeed predicts a production rate within this range and LFV should be observed. As 
discussed above, the background for such lepton-flavor violating processes is relatively clean 
and the signal-to-noise is high, so that observation would be more optimistic. 

As we know, to judge validity of a new physics model, usually one cannot decide its 
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existence only by a unique experiment, but needs a combination of data from various mea- 
surements. Such a synthesis analysis may help to confirm a new model. Phenomenological 
analysis on effects induced by unparticle physics has been carried out for various production 
and decay processes and even the cosmology. In general, the concerned parameters, such as 
its coupling to SM matter, scale dimension dy and other characteristics are constrained by 
earlier analysis, therefore it seems to be the time to let experiment tell us the validity of the 
whole theory. Thus our conclusion is that observation of the lepton-flavor violating decays 
may provide an important test of the unparticle scenario. If the observation confirms the 
predicted value, the unparticle theory would be greatly supported (even though relatively 
large uncertainty still remains), by contrary, if the reaction is not seen at BESIII, one may set 
a more stringent constraint on the parameters which would be applied to other experiments 
for further tests. 
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